Anisotropic Hydrodynamics for Rapidly Expanding Systems 



o 

(N 

Oh 
< 



i 



> 

in 

O 
O 



X 

c3 



Wojciech Florkowski, 1,2 Radoslaw Ryblewski, 2 and Michael Strickland 3,4 

1 Institute of Physics, Jan Kochanowski University, PL-25J h 06 Kielce, Poland 

2 The H. Niewodniczanski Institute of Nuclear Physics, 

Polish Academy of Sciences, PL-31342 Krakow, Poland 

3 Department of Physics, Kent State University, Kent, OH 44^42 United States 

4 Frankfurt Institute for Advanced Studies, Ruth-Moufang-Strasse 1, D-60438, Frankfurt am Mam, Germany 

(Dated: April 3, 2013) 

We exactly solve the relaxation-time approximation Boltzmann equation for a system which is 
transversely homogeneous and undergoing boost-invariant longitudinal expansion. We compare the 
resulting exact numerical solution with approximate solutions available in the anisotropic hydro- 
dynamics and second order viscous hydrodynamics frameworks. In all cases studied, we find that 
the anisotropic hydrodynamics framework is a better approximation to the exact solution than 
traditional viscous hydrodynamical approaches. 



The application of relativistic viscous hydrodynamics 
is important in a wide variety of situations including, 
for example, the dynamics of high energy astrophysical 
plasmas and the quark-gluon plasma created in relativis- 
tic heavy-ion collisions. Since the seminal work of Israel 
and Stewart [TJ [2] there have been many papers that 
have addressed the questions of how to apply and sys- 
tematically improve relativistic viscous hydrodynamics 
[3ti2Tj . Recently, a new framework called anisotropic hy- 
drodynamics (aHydro) has emerged for describing the 
non-equilibrium dynamics of relativistic systems [2"2H2"5] . 

In contrast to conventional viscous hydrodynami- 
cal treatments, aHydro does not implicitly rely on 
an assumption that the system is approximately 
isotropic in momentum-space. Instead, momentum-space 
anisotropies are built in at leading order by utilizing a 
spheroidal form for the one-particle distribution function. 
For systems that are boost invariant and homogeneous 
in the transverse direction it has been proven that the 
aHydro framework reduces to second order viscous hy- 
drodynamics in the limit of small anisotropies |23j . In 
addition, the framework reproduces the longitudinal free- 
streaming limit and in all cases the one-particle distribu- 
tion and the transverse/longitudinal pressures are guar- 
anteed to be positive. 

Given these appealing features, the aHydro framework 
seems to offer a promising alternative to standard viscous 
hydrodynamical treatments. However, in order to judge 
the efficacy of the framework it is desirable to have an 
exactly solvable case with which to compare the various 
approximations. With this in mind, in this Letter we 
exactly solve the Boltzmann equation for a transversely 
homogeneous boost-invariant system of massless particles 
in the relaxation time approximation. We then compare 
the resulting exact solutions with aHydro and two differ- 
ent second order viscous hydrodynamics approximations. 

Our starting point is the Boltzmann equation 
P >1 9 fl f(x,p) = C[f(x,p)} where f(x,p) is the one-particle 
distribution function and C[f] is the collisional kernel. 
We will focus on the case of the relaxation time approx- 



imation (RTA) 



C[f} = 



P^ 



• eq 



f c ^\P^,T{xyj - f(x,p) 



(1) 



where u^ is the local rest frame four velocity, r oq is the 
relaxation time which may depend on proper time, and 
/ C q is an equilibrium distribution function that may be 
taken to be a Bose-Einstein, Fermi-Dirac, or Boltzmann 
distribution. The effective temperature T{t) appearing 
in the argument of the equilibrium distribution function 
is fixed by requiring dynamical energy-momentum con- 
servation [5Tj . 

We focus on boost invariant systems which are ho- 
mogenous in the transverse direction in which case the 
dynamical variables only depend on the proper time. We 
additionally specialize to the case of massless particles. 
Defining variables w = tpr, — zE and v = Et — p^z |32l - 
134] . where pi, is the particle momentum along the z di- 
rection, the left hand side of the Boltzmann equation can 
be written simply as p^d^f = (v/r)d T f. This allows one 
to solve the RTA Boltzmann equation exactly 

f(T,w,p x ) = D(T,T )f (w,p±) 

-^—D(r,r')f cq (r',w, Pl _),(2) 

where To is the initial proper time, /o is the ini- 
tial non-equilibrium distribution function, and 



is the damping 



C(r 2 ,r 1 )=exp[-/; 2 dr"r- 1 (r")_ 
function. This solution is similar to the one obtained 
originally by Baym [31] , see also [35l [36] . We have gen- 
eralized it to an arbitrary initial condition at r 7^ and 
allowed for the possibility that the equilibration time r eq 
is time dependent. In the relaxation time approximation 
one finds r cq = 5r]/(TS) where rj is the shear viscosity, S 
is the entropy density, and T is the effective temperature 
which we will specify below [37l [38]. Herein we will 
assume that rj/S is time independent. 
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FIG. 1. (Color online) Pressure anisotropy as a function of 
proper time assuming ft = and To = 600 MeV at to — 0.25 
fm/c for Anrj/S — 1 (top), 3 (middle), and 10 (bottom). Exact 
solution (black solid), aHydro (AH) approximation (red long- 
dashed), Israel-Stewart (IS) approximation (blue dot-dashed), 
and full second order (DNMR) approximation (brown dotted) 
|20| are compared. 



Based on Eq. ^, one can evaluate the energy density 
via £ (t) = g J dPv 2 f(T,w,p±)/T 2 where g is the degen- 
eracy factor and dP = 2d A p5(p 2 )8(p a ) = v^ 1 dwd 2 pT- 
Integrating Eq. pi one obtains an integral equation for 
the energy density 



£{r) = D(t,t )1I(Z fs (t))/K(Z ) 

"'' D(r,r')£(r')n((^) 2 



t 'eq 



(T>) 



1 , (3) 



where £ = £/£q is the energy density scaled by the 
initial energy density, ft) is the initial momentum- 
space anisotropy, £fs( t ) = (1 + (,o)( T / T o) 2 — 1 ; an d 
K(z) = \ [(1 + z)- 1 + arctan(v5)/ 1 /i] ■ 

Equation pi can be solved numerically using the itera- 
tion method. From the resulting energy density, one can 
solve for the effective temperature via £{t) — jT 4 (t) 
where 7 is a constant which depends on the particular 
equilibrium distribution function assumed and the num- 
ber of degrees of freedom. The resulting effective temper- 



FIG. 2. (Color online) Pressure anisotropy as a function of 
proper time assuming ft = 10 and To = 600 MeV at tq = 
0.25 fm/c for Airq/S = 1 (top), 3 (middle), and 10 (bottom). 
Labeling is the same as in Fig. [l] 



ature allows one to determine the distribution function 
/ eq at all proper times and, with this, the full particle 
distribution function can be obtained using Eq. pi . Ad- 
ditionally, one can determine the number density, longi- 
tudinal pressure, and transverse pressure by integrating 
the distribution function multiplied by v/t, w 2 /t 2 , and 
Pt/2, respectively. 

We will compare the resulting numerical solutions with 
the solution of the aHydro (AH) equations which are 
obtained by evaluating the zeroth and first moments 
of Boltzmann equation in the relaxation time approx- 
imation assuming a spheroidal form for the distribu- 
tion function /ah = /eq(((P/xU M ) 2 + £(JV^) 2 )/A 2 J with 
u» = (t,0,0,z)/r and z^ = (z,0,0,t)/r. With this as- 



sumption, one obtains [23 



1 



1 T C 2 



A T. 
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AH 



i-^ 3/4 (oyrTe 
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where A 



d T £ + T d T A 
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n AH 



the transverse temperature 
5rj/(2AS) is the relaxation time |39j . 



.(4) 

and 
We 
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FIG. 3. (Color online) Pressure anisotropy as a function of 
proper time assuming £o = and To = 300 MeV at to = 
0.25 fm/c for 4-irrj/S = 1 (top), 3 (middle), and 10 (bottom). 
Labeling is the same as in Fig. IT] 



FIG. 4. (Color online) Pressure anisotropy as a function of 
proper time assuming £o = 10 and To = 300 MeV at tq = 
0.25 fm/c for 4Tvrj/S = 1 (top), 3 (middle), and 10 (bottom). 
Labeling is the same as in Fig. [TJ 



note that in the limit £ —¥ one has A — > T and 
t^ h = r cq /2. The time evolution of £ and A is obtained 
by solving Eqs. Q and using these one can straight- 
forwardly compute the time dependence of the energy 
density, transverse pressure, longitudinal pressure, and 
number density [53J HO] . 

In addition we will compare the exact solution with two 
second order viscous hydro prescriptions, both of which 
can be written compactly as 



d T £ 



d T U 



£ + V 



n 



r r 

n An 
— +,— 







n 



(5) 



where II is the shear and t„. = 5n/(TS) is the shear relax- 
ation time. In the majority of the literature, practitioners 
use = 4/3 which we will refer to as the Israel-Stewart 
(IS) prescription. We will also compare the exact so- 
lutions with the complete second order treatment from 
Ref. |20j which, within the relaxation time approxima- 
tion, gives — 38/21. We will refer to the second choice 
as the DNMR prescription 41j. In both cases one can 



compute the transverse pressure via Vt = V + LT/2 and 
the longitudinal pressure via Vl = V — II. To be consis- 
tent with the exact solution and the aHydro approxima- 
tion we assume an ideal equation of state for the visous 
hydrodynamical approximations. 

We now turn to our results. For all results shown we 
have assumed that the initial distribution function was 
spheroidal in form but for the exact solution we do not 
restrict the distribution function in any way after this 
point in time. In Fig.fTJwe show the pressure anisotropy 
as a function of proper time assuming an initial isotropic 
plasma with £ = and T = 600 MeV at r = 0.25 
fm/c. The three panels show three different assumed val- 
ues of the shear viscosity to entropy ratio corresponding 
to Attti/S <E {1,3,10}. In the figure the aHydro, Israel- 
Stewart, and full second order viscous hydrodynamics ap- 
proximations |20] are compared with numerical solution 
of the Boltzmann equation. In Fig. [2] we plot the same 
quantity for a system possessing an initial momentum- 
space anisotropy corresponding to £ — 10- In Figs. [3] 
and [4] we present the pressure anisotropy subject to the 
same initial conditions and values of 77/ S only changing 
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FIG. 5. (Color online) Distribution function contour lines re- 
sulting from the exact numerical solution (black solid) and the 
aHydro approximation (red dashed) at three different times 
r = 0.25 fm/c (top), 1 fm/c (middle), and 10 fm/c. The pa- 
rameters used where £o = 10 and To = 600 MeV at to = 0.25 
fm/c for 4-ktj/S = 3. 



the initial (effective) temperature to To = 300 MeV. 

As these figures demonstrate, the aHydro approxima- 
tion is always closer to the exact solution than the IS 
and DNMR approximations. The IS approximation is 
the worst approximation to the exact solution in all cases 
shown and for the case Atttj/S — 10 it even predicts a neg- 
ative longitudinal pressure for the majority of the time 
shown. The DNMR approximation represents a signif- 
icant improvement over the IS approximation; however, 
we note that if one increases the shear viscosity to entropy 
ratio even further, the DNMR approximation also pre- 
dicts negative longitudinal pressures. Within the aHy- 
dro approximation, on the other hand, the pressures are 
guaranteed to be positive at all times. In addition, within 
aHydro the one-particle distribution function is guaran- 
teed to be positive at all times. 

In Fig. [5] we show three contour plots which compare 
the evolution of the one-particle distribution function 
obtained from the numerical solution of the Boltzmann 
equation with the aHydro approximation. As this fig- 
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FIG. 6. (Color online) Particle production measure A„ = 
r//ron(r/)/n(ro) — 1 as a function of 4irr]/S. Exact solu- 
tion (black points), aHydro (AH) approximation (red long- 
dashed), Israel-Stewart (IS) approximation (blue dot-dashed), 
and full second order (DNMR) approximation (brown dotted) 
[20] are compared. 



ure shows, the exact solution shows deviations from the 
spheroidal form of aHydro at early times. The magnitude 
of these deviations decreases as the shear viscosity to en- 
tropy ratio decreases. However, as Figs. l]4 demonstrate, 
moments of the distribution function are described quite 
well. We note that the shape of the contours shown in 
Fig. [5] for the exact result suggest that it may be more 
effective to describe the distribution function in terms of 
a linear superposition of two spheroidal forms, one which 
is related to the free streaming part contribution from 
the initial distribution function and the other which de- 
scribes the late-time dynamics. 

In Fig. [6] we show the particle number generation mea- 
sured via A n = t//to n(r/)/n(ro) — 1 where 17 is the 
freeze-out temperature defined by when the effective tem- 
perature drops below Tf = 150 MeV. In the figure, we 
compare A n obtained via numerical solution of the Boltz- 
mann equation, the aHydro approximation, the IS vis- 
cous hydrodynamical equations, and the DNMR viscous 
hydrodynamical equations. On physical grounds one ex- 
pects this ratio to vanish in the ideal hydrodynamical 
limit (ji/S —¥ 0) and the free-streaming limit (r//S — > 00). 
As this figure demonstrates both the IS and DNMR ap- 
proximations predict that A n is a monotonically increas- 
ing function of the shear viscosity to entropy ratio. On 
the other hand, although the aHydro framework has too 
much particle production at large rj/S, it has the right 
qualitative behavior. 

In conclusion, we have presented an exactly solvable 
case in which the Boltzmann equation can be straight- 



forwardly solved numerically. This exact solution allows 
one to test the accuracy of different dynamical approxi- 
mation schemes for a variety of initial conditions, values 
of the shear viscosity, etc. We demonstrated that the 
aHydro approximation is closer to the exact solution than 
both the IS and DNMR approximations in all cases. This 
is remarkable since the aHydro equations themselves are 
essentially zeroth order in a general anisotropic expan- 
sion of the one-particle distribution function, whereas the 
other approximations are at their respective second order 
of approximation. If deviations from the exact spheroidal 
form are taken into account, one expects the quality of 
the aHydro approximation to further improve. We note 
in closing that the exact solutions obtained herein can 
perhaps be used to create a better approximation frame- 
work which takes into account such deviations in a sys- 
tematic manner. 
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